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Abstract. We present a necessary and sufficient condition for a 3 X 3 matrix 
to be unitarily equivalent to a symmetric matrix witli complex entries, and 
an algoritlim whereby an arbitrary 3x3 matrix can be tested. Tliis test 
generalizes to a necessary and sufficient condition tfiat applies to almost every 
n X n matrix. Tlie test is constructive in that it explicitly exhibits the unitary 
equivalence to a complex symmetric matrix. 



1. Introduction 

Definition. A complex symmetric matrix is an n x n matrix T with complex entries 
such that T — . We wiU often refer to such a matrix as a CSM. We will often 
refer to a matrix that is unitarily equivalent to a complex symmetric matrix as 
being UECSM. 

Definition. A conjugation is a conjugate-linear operator C : C" ^ C" which is 
both involutive (C^ = /) and isometric {{x,y) — {Cy,Cx) for all x,y £ C"). We 
say that a linear operator T on C" is C-symmetric if T = CT*C. 

The matrices that are unitarily equivalent to complex symmetric matrices can 
be characterized in terms of C-symmetry. We recall a consequence of [B], Prop. 2]: 

Proposition 1. An nx n matrix T is unitarily equivalent to a complex symmetric 
matrix if and only if there exists a conjugation C for which the linear operator 
induced by T is C-symmetric. 

For more on complex symmetric operators, see [HIT]. It is already known that 
all 2 X 2 matrices are unitarily equivalent to a CSM ([SI Ex. 6], [2 Cor 3.3]) but 
that not all 3 X 3 matrices are p", Ex. 7]. It is also known that all matrices are 
similar to a CSM (j6j Ex. 4], [HI Thm. 4.4.9]), which makes it difficult to determine 
if a particular matrix is unitarily equivalent to a complex symmetric matrix. For 
instance, exactly one of the following matrices is unitarily equivalent to a CSM: 



Ti 







To 







For more, see Example [H 



This note presents necessary and sufficient conditions for a 3 x 3 matrix to be 
UECSM, as well as an algorithm with which one can apply these conditions to an 
arbitrary 3x3 matrix. This test also applies to almost every n x n matrix T, and 
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if T is UECSM then it possible to explicitly construct a conjugation C for which T 
is C-symmetric. 

2. Eigenvectors of the Cartesian decomposition 

Any matrix T has a unique Cartesian decomposition T = A + iB, where A and 
B are self-adjoint. A and B are given explicitly by the formulas A= i(T + T*) 
and B = -^{T — T*). It follows that T is C-symmetric for some conjugation C if 
and only if A and B are C-symmetric. 

We recall a result from [6l Lemma 1]: 

Lemma 1. IfC is a conjugation on C", then there is an orthonormal basis {e^}"^]^ 
such that Cci — e,; for all i. 

Such an orthonormal basis is called C-real. The following lemma shows that 
one way such bases arise is from orthonormal bases of eigenvectors of self-adjoint 
matrices. 

Lemma 2. If A is a C-symmetric self-adjoint matrix, then there exists a C-real 
orthonormal basis of eigenvectors of A. 

Proof. If A is an eigenvalue of yl with eigenvector x, then ACx — CAx = CXx = XCx, 
since A is real. Since C is involutive, the above equation yields that C must map 
the eigenspaces of A onto themselves. Thus we can decompose C = Ci ® • • • ® Cm 
where the Ci are conjugations on the eigenspaces of A corresponding to distinct 
eigenvalues. Since the ith eigenspace must have a C^-real orthonormal basis, and 
the eigenspaces of A are mutually orthogonal, we are done by Lemma [TJ □ 

The following lemma characterizes all UECSM matrices in terms of such or- 
thonormal bases of eigenvectors. 

Lemma 3. If T = A -\- iB is an n x n matrix, then T is unitarily equivalent to a 
symmetric matrix if and only if there exist orthonormal bases of eigenvectors {ci} 
and {fi} of A and B respectively such that {ei,fj) e M for all 1 < i, j < n. 

Proof. First assume that T is UECSM. By Proposition [1] there must exist a con- 
jugation C on C" such that CT*C = T, and thus CAC = A and CBC = B. By 
Lemma [21 A and B have C-real orthonormal bases of eigenvectors {ci} and {fi}, 
respectively. For these we have, 

(e»,/j) = {Cfj,Cei) = (fj,ei) = {e,Jj), 

and so (e^, fj) £ R. 

Conversely suppose there exist such {ei} and {fi}. Then define C : C" ^ C" by 

n 

Cx = ^ (x, ei)ei. 

i=l 

It is easy to check that C is a conjugation, and we can calculate 

n n 
Ccj — ^ ^ (^J 7 — ^ ^ ^ijj^i — 
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and, 



i=\ i=l 

Since we know Aci = XiCi for some A; € R, we can calculate, 

CACei = CAci = XiCi = Aci. 

By linearity, it follows that Ax = CACx for all x e C". Similarly one can show 
that CBC = B, and so T is C-symmetric. □ 

Remark. In the previous proof, we could have let C be complex conjugation with 
respect to {/j} instead. 

Although the following corollary is well-known [4, Ex. 2.8], the above lemma 
provides another proof. 

Corollary 1. Every normal matrix is unitarily equivalent to a complex symmetric 
matrix. 

Proof. T = A + iB being normal implies that A and B share a basis of eigenvectors, 
and so in the hypothesis of Lemma [3] we may pick {ci} and {fi} to be the same 
orthonormal basis. □ 

Corollary 2. Let T = A + iB be an n x n matrix that is unitarily equivalent to a 
complex symmetric matrix, and let {ci} and {fi} be a pair of orthonormal bases of 
A and B, respectively, as in Lemma\^ Then for any A' and B' such that the {ei} 
are eigenvectors for one and {fi} are eigenvectors for the other, T' — A' + iB' is 
unitarily equivalent to a complex symmetric matrix. 

3. A TEST FOR n X n matrices 

Definition. We say that a pair of orthogonal bases {gi} and {/ij} for C" are proper 
if (gi, /ii) e M and hj) = =^ 1^1 and j ^ 1. If we let M = {{g,, 
then this is equivalent to the top-left entry of M being real while the first row and 
column contain no zeros. 

Theorem 1. Let T = A + iB be an n x n matrix, and let {gi} and {hi} be any 
proper pair of orthogonal bases of eigenvectors of A and B, respectively. Then T is 
unitarily equivalent to a complex symmetric matrix if for 2 < i, j < n, 

. e M. (1) 

{9z,hi){gi,hj} 

Moreover, if A and B both have n distinct eigenvalues, then ^ is also necessary 
for T to be UECSM. 

Proof. Suppose that llj holds. Define ei = gi, f\ — h\ and otherwise. 



Once normalized, these bases satisfy Lemma [H showing that T is UECSM. 

Conversely, suppose that T is UECSM and that A and B have n distinct eigen- 
values. By Lemma [HI there are orthonormal bases of eigenvectors {ei} and {fi} of 
A and B, respectively, such that {a, fj) G M. As the eigenspaces of A and B are 



4 



JAMES E. TENER 



one-dimensional, we can reorder these bases so that gi — uJiCi and hj ~ Qfj for 
unimodular uji, (j G C. Then for 2 < i, j < n, 



which is real since each {ei, fj) G M. The pair {gi} and {hi} being proper ensures 
that {gi,hi) G M \ {0}, and so the preceding equation yields that {gi} and {hi} 



Remark. The condition ([T]) of Theorem [T] can be visualized using matrices. If {ei} 
and {fi} are a proper pair of orthogonal bases of eigenvectors, we can consider the 
matrix M = (rriij) = ((ci, fj)), which will be unitary if both bases are normalized. 
Thinking of M with 1 x (n — 1) blocking. 



M = 







Cl 


D 



the condition says that each element in the lower-right block D has the same ar- 
gument as the product of the first element in its row and the first element in its 
column. 

Remark. If a matrix is verified to be UECSM via Theorem [U then an explicit 
conjugation C is provided in terms of the eigenvectors of A and B via the proofs of 
Lemma[3]and Theorem[TJ Let {ci} and {fi} be the orthonormal bases from Lemma 
|3l [/ be the transition matrix from the standard basis to either {ci} or {fi}, and 
let J be the coordinate- wise complex conjugation operator on C". It follows from 
[i Sec. 3.2] that U*TU is a CSM and that 

C = UJU* = UU^J. 

is a conjugation with respect to which T is C-symmetric. 

While there are already multiple proofs that all 2 x 2 matrices are UECSM ([6l 
Ex. 6],nn, Cor S.S],!?", Cor. 1]), we can use Theorem [1] to provide another proof 
while explicitly exhibiting a conjugation C for which T = CT*C . 

Corollary 3. Every 2x2 matrix is unitarily equivalent to a complex symmetric 
matrix. 

Proof. Let T = yl + zi3bea2x2 matrix. We may assume without loss of generality 
that A and B have 2 distinct eigenvalues, and that they do not share an eigenvec- 
tor, as in either case T would be normal and therefore UECSM by Corollary [1] 
Let {e^} and {fi} be orthonormal bases of eigenvectors of A and B, respectively. 
Compute the unitary U — {ui_j) — {{ei,fj)), and observe that since A and B do 
not share an eigenvector, no entry of U is equal to 0. Multiply ei by a unimod- 
ular constant so that (ei,/i) G M. Since the columns of U must be orthogonal, 
iti_iUi.2 + U2,iU2,2 = which yields 



{9i,hj){gi,hi) 
{9i,hi){gi,hj) 



(^»e,,Cj7j)(^iei,Ci/i) 
{uJie,,Cifi){ujiei,Qfj) 

^»^iCjCi(e»:/?)(ei:/i) 
^i^^iCjCi{ei,fi){ei,fj) 
(ej,/j)(ei,/i) 



(e»,/i)(ei,/j) 



satisfy ([T]). 



□ 



W2,2 



Ul,2W2,l 




UNITARY EQUIVALENCE TO A COMPLEX SYMMETRIC MATRIX 



5 



As ui 1 e R, it follows that e M. By Theorem [H T is UECSM. If we 

multiply 62 by a unimodular constant so that (e2, /i) G M and let y = [ ei | 62 ] , 
then T is C-symmetric for Cx — VV*x by the remark following Theorem [1] □ 

4. An ALGORITHM FOR 3x3 MATRICES 

In this section we introduce an algorithm for determining whether or not a 
given 3x3 matrix is unitarily equivalent to a complex symmetric matrix. We first 
require a few preparatory results. The following proposition allows us to easily 
answer afiirmatively in certain cases. 

Proposition 2. Let T = A + iB be a 3 x 3 matrix. If either of the following 
conditions hold, then T is unitarily equivalent to a complex symmetric matrix: 

(i) A or B has a repeated eigenvalue 

(ii) A and B share an eigenvector 

Proof First (i). Note that T is UECSM if and only if A/ is. Ii B has a repeated 
eigenvalue A, then B — XI either has rank or 1. If B — XI — 0, then T — iXI = A, 
which is self-adjoint and therefore UECSM by Corollary [1] li B — XI has rank 1, we 
can show that T is UECSM by simplifying an earlier result 8, Cor. 5] for the finite 
dimensional case. Let w be a unit vector that spans Kan{B — XI), and let Ai be the 
cyclic subspace generated by v under A. Clearly M is invariant under A, and it is 
easy to check that A4-^ is as well. T \j^± is self-adjoint and T \^ is cyclic self- 
adjoint, so by the Spectral Theorem we can assume without loss of generality that 
V = (1, 1, 1) and that A is diagonal [21 Thm. 2.11.2]. Since Ran(B- A/) = span{i;}, 
we have 

B = b 

for some 6 e K. This yields that T = A + iB is UECSM. The proof is similar if A 
has a repeated eigenvalue. 

Now (ii). If A and B share an eigenvector, then up to unitary equivalence we 
know that A = Ai (B A2 and B = Bi® B2 where Ai and Bi are 1x1 matrices and 
A2 and B2 are 2 X 2 matrices. This yields that T = (Ai + iBi) {A2 + 1-62), and 
because aU 1 x 1 and 2 x 2 matrices are UECSM, T is UECSM. □ 

The next lemma tells us that when Proposition[2]does not imply that T = A + iB 
is unitarily equivalent to a CSM, it is easy to construct a proper pair of orthogonal 
bases to which one can apply Theorem [TJ 

Lemma 4. If T = A + iB is a 3 x 3 matrix which does not satisfy either hypothesis 
of Proposition\^ then any two orthogonal bases of eigenvectors {ci} and {fi} of A 
and B, respectively, can be made proper by reordering them and scaling ei. 

Proof. Without loss of generality we may assume that ||ei|| = ||/i|| = 1, since the 
conditions of being proper are not affected by multiplying the basis vectors by 
real scalars. Note that no is orthogonal to more than one fj, since otherwise 
it would be a scalar multiple of the third element of {fj}, a contradiction. Simi- 
larly no fixed is orthogonal to more than one e^. In terms of the unitary matrix 



111 
111 
111 



6 



JAMES E. TENER 



U — (uij) — {{ci, fj)), this means that no row or column has more than one 0. 

We claim that there is at most one in U. If U had more than one 0, they must 
be in different rows and columns, so we could reorder the bases so that 



U = 



* * 
* * 
a 6 * 



where the *'s represent arbitrary complex numbers. To preserve the orthogonality 
of columns, we must have a = or 6 = 0. In either case, there must be more than 
one in a single column, which we have already excluded. By reordering the bases, 
we can ensure that the entry is not in the first row or column. If (ei, /i) ^ R, 
multiply ei by and then {ei} and {fi} will be a proper pair. □ 

Using the preceding results, we can construct an algorithm that will decide 
whether or not a 3 x 3 matrix T is unitarily equivalent to a complex symmetric 
matrix. Since none of the operations are more complicated than finding roots of 
cubic polynomials, it can be performed using exact values, assuming the data is 
given exactly. The author implemented it in Mathematica without much difficulty. 
The algorithm is: 



Algorithm. Given a 3 x 3 matrix T, 



(1) Compute A = i(T + T*) and B = j-{T - T*). 

(2) Compute the eigenvalues of A and B. If either A oi B has a repeated 
eigenvalue, then T is UECSM by Proposition [2l 

(3) Compute arbitrary sets of eigenvectors {gi} and {hi} of A and B, respec- 
tively, and compute the matrix M = {mi_j) = {{gi, hj))ij. 

(4) If M has more than one entry equal to 0, then T is UECSM (Lemma 21 
Proposition [5]) . Otherwise, reorder the rows and columns of M so that the 
entry is not in the first row or column (so that {g,} and {hi} form a 
proper pair). Scale gi by 

(5) By Theoremlll T is UECSM if and only if for all 2<i,j <3, 



(6) If T is UECSM, one can exhibit a corresponding conjugation C by first 
normalizing {gi} and scaling 172 and g^ so that {gi,hi) € M for all i. If 
= [ 51 I 52 I 53 ] , then by the remark following Theorem [1] U*TU is 
complex symmetric and T is C-symmetric with respect to the conjugation 
Cx = UU*x. 

It is worth noting that steps 1, 3 and 5 carry through to the nxn case as long as 
A and B have n distinct eigenvalues and a proper pair of bases can be found. Step 
4 is no longer valid, as for n > 3 the preceding conditions do not guarantee that T 
is UECSM. A generalization of this algorithm to n x n matrices will be discussed 
in Section [5] The following examples illustrate the steps of the algorithm. 
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Example 1. Theorem [T] provides another proof of the fact from [8i Ex. 1] that for 
a, 6 7^ the matrix 

"060 
r= a 


is unitarily equivalent to a complex symmetric matrix if and only if \a\ = By 
dividing by 6, it is enough to consider matrices of the form 

"010 

r= a 


For this T, one can verify that 

1 



A 



5 

-2 f 

I 



and B 



— 

" 2 



Furthermore, the eigenvalues of both A and B arc {0, + |ap, — i-y^l + |ap}, 
the eigenvectors of A are 





a 






1 




1 


91 = 







, 92 ^ 






, 53 = 






-1 






a 




a 


eigenvectors of B are 










a 






1 




1 


hi = 

















1 






—a 




—a 



The matrix M — {{gi, hj))f required by the algorithm is given by: 



M = 



l-|a|2 2 2 
2 ^ ^ 
2 (3 



where f3 — 1 + If |a| 7^ 1, then the bases {gi} and {hi} are proper and since 

/3 has non-zero imaginary component, T is not UECSM by Theorem [TJ 



If \a\ = 1, then (3 — 2i. We can relabel the vectors of each basis and scale the 
new ei by —i to get the matrix, 

"1 -1 

M' = 2 -i i 1 
1 1 

e M for i, j > 2, so T is UECSM by Theorem [T] 



It is easy to check that — — 
Example 2. Let 

' l + ii {-2-i)V2 -l-4i 
r= ^^/2 iV2 

-1 {2-i)V2 1 

In this example, we prove that T = A + iB is unitarily equivalent to a complex 
symmetric matrix and use the full algorithm to find a conjugation C with respect 
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to which T is C-symmetric. Per step 2, we first calculate the eigenvalues of A and 
_B, which are 

{2(1 + V2) ,-2,2(1- V2)} and |2(l + V3) ,2(1- Vs) ,0}, 

respectively. Neither A nor B has a repeated eigenvalue, and so we calculate the 
eigenvectors of A: 



91 

and the eigenvectors of B: 

hi = 



' -1 - 2iV2 ' 




1 




" -1 + 2iV2 ' 


2 + iV2 


, 92 = 






-2 + iV2 


3 




1 




3 







r -1 + — 1 

^/3 




1 




, h2 = 


1 


, = 


i^/2 


1 






1 



Since A and B do not share an eigenvector, we must use the full algorithm and 
not a shortcut provided by Proposition [2l Next we compute AI ~ {{gi,hj)): 



M 



{2 + iV2) {3 + VS) -fi(-2i + V2) (-3 + %/3) 4-4iV2 



4 

V3 



(2 - iV2) (3 + V3) 



4 



(2i + V2) (-3 + V3) 4 + 4i\/2 



If we let a = ^-^^ and /3 
{hi, h2, -ih^} we get. 



l™3.l| 
ITI3 



^, rewriting M with respect to {agi, g2, ^93} and 



2{V2 + VQ) 2\/2(-1 + V3) 4\/3 

4_ J_ Q 

\/3 V3 " 

2(\/2 + \/6) 2\/2(-l + V3) -4\/3 



As all of the entries are real. Theorem 1 says that T is UECSM. Letting {ei} 



{||ff7|| , nflii , nlffli } and U = [eilcaleg], we can construct a conjugation C for which 
T is C-symmetric and a complex symmetric matrix that T is unitarily equivalent 
to: 





1 


i 


1 " 








2 




2 





















2 


z 

V2 


2 







U*TU 



Example 3. Let 





' l + V2 + i -i 


i 


= 2 




-1 


—i 




i 


—i 


1 -V2 + i _ 


" 


7 




"0 7 " 





1 -5 


, T2 = 


1 -5 





6 




3 
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Applying the same method as in the previous example yields that Ti is C-symmetric 
where 





Xi 




c 








X3 





6(-19+6iy74) 

3(J25 
42(l9-6jV74) 

3025 

fe(l9-6*%/74) 



42(l9-6i\/74) 

3025 
19(-19+6iV74) 



19 - 6iV74 



605 



6 
605 



3025 
(19 - 6iVU) 



7 

605 

^^(19-6*^74) 
6§5(-19 + 6*V74) 



and that Ti is unitarily equivalent to 



56 
37 



55 35V55 
'37 74 



55 
37 
35\/55 
74 



56 
37 



+ i 



35V55 
74 



37 35V55 
74 
147 
37 



We also get that T2 is not UECSM. It is easy to check that if T2 has Cartesian 
decomposition T2 — A + iB, both A and B have 3 distinct eigenvalues and that they 
do not share an eigenvector. None of the eigenspaces of A are orthogonal to any 
of the eigenspaces of B, and so it is easy to construct a proper pair of orthogonal 
bases of eigenvectors of A and B. It is also easy to check that these will not satisfy 
condition ^ of Theorem [1] 



5. Applications for generic n x n matrices 

Theorem [T] provides a sufhcient condition for an arbitrary n x n matrix to be 
unitarily equivalent to a complex symmetric matrix, so long as a proper pair of 
bases can be found, and if A and B each have n distinct eigenvalues then Theorem 
[T]is necessary as well as sufficient. It is not guaranteed that an arbitrary pair of or- 
thogonal bases can be made proper simply by reordering and scaling the elements. 
However, for a random matrix T, its Cartesian components A and B will both 
have n distinct eigenvalues with probability one [TOl Chapter 3]. Moreover, with 
probability one an arbitrary choice of bases of eigenvectors of A and B can easily 
be made proper, as in Lemma SI For more on random self-adjoint matrices, one 
can consult [lOj . In light of the above, the algorithm given in the previous section 
will almost surely apply to a random nx n matrix, which is useful for probabilistic 
searches. 

For instance, in [51 Sec. 5] the authors determine that any 4x4 partial isometry 
T with RankT = 1 or RankT = 3 is UECSM. T is also trivially UECSM in the 
cases where RankT = or RankT = 4. However, they were unable to answer 
whether or not every 4x4 partial isometry with Rank T = 2 is unitarily equivalent 
to a complex symmetric matrix. All such partial isometrics are unitarily equivalent 
to a matrix of the form UP, where P the projection onto the first two standard 
basis vectors and U is unitary. We used Mathematica to test 100,000 matrices of 
the form U T, with random unitary components generated by taking the matrix 
exponential of a random skew-Hermitian matrix. All of them were UECSM, and 
so we conjecture: 

Conjecture 1. Every rank-two 4x4 partial isometry is unitarily equivalent to a 
complex symmetric matrix. 
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